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Absence of Anomalous Tunneling of Bogoliubov Excitations for Arbitrary Potential 
Barrier under the Critical Condensate Current 
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We derive the exact solution of low energy limit of Bogoliubov equations for excitations of 
Bose-Einstein condensate in the presence of arbitrary potential barrier and maximum current 
of condensate. Using this solution, we give the explicit expression for the transmission coeffi- 
cient against the potential barrier, which shows partial transmission in the low energy limit. 
The wavefunctions of excitations in the low energy limit do not coincide with that of the con- 
densate. The absence of the perfect transmission in the critical current state originates from 
local enhancement of density fluctuations around the potential barrier. 

KEYWORDS: Bose-Einstein condensate, Gross-Pitaevskii equation, Bogoliubov equations, anomalous tun- 
neling, critical current 



Anomalous tunneling, which was first discovered 
by Kovrizhin and his collaborators,^-' is an interesting 
theoretical prediction concerning the tunneling proper- 
ties of Bogoliubov excitations of weakly-interacting Bose 
gas in superfluid phase; Bogoliubov excitations tunnel 
across a potential barrier without reflection in low energy 
limit. This prediction has a relevance to experiments in 
magnetically trapped dilute Bose-Einstein condensates, 
where the propagation of collective excitations^-* and the 
structure factor of density fluctuation'^' have been ac- 
counted for by Bogoliubov theory.^^ 

Generalized or related issues of the anomalous tun- 
neling have been considered by several authors. Par- 
ticularly, Danshita et al.^^ have considered the tunneling 
properties of excitations of Bose system in the presence 
of supercurrent of condensate, and have found that (a) 
the perfect transmission occurs even when the conden- 
sate current exists, except for the critical current state; 
(b)under the critical current, the perfect transmission 
disappears, and only a partial transmission occurs. 

Recently, the anomalous tunneling has been proved 
for an arbitrary barrier in the absence of condensate cur- 
rent, using the fact that wavefunctions of Bogoliubov 
excitations in the low energy limit coincide with that of 
a condensate. Subsequently, Ohashi-Tsuchiya discussed 
the origin of anomalous tunneling considering the con- 
densate current'^^) and attributed (a) to the similarity 
between a low-energy Bogoliubov excitation and a con- 
densate wavefunction. 

From the scenario, the anomalous tunneling is 
expected to be a universal phenomenon related to excita- 
tions of superfluid condensate. Therefore the anomalous 
tunneling should be understood within the framework 
applicable to arbitrary shape of potential barrier. How- 
ever, earlier works^'^'"'^^-' but refs. 10, 11 have been based 
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on the exact solutions for a particular shape of poten- 
tial barrier and there is no discussion on (b) for general 
potential barrier so far. It is thus important to clarify 
whether the scenario^ ^ ■ "'^^^ of the origin of the anomalous 
tunneling explains consistently both the perfect trans- 
mission in the non-critical states and the partial trans- 
mission in the presence of critical condensate current for 
general potential barrier. 

In this Letter, we derive the exact expression for 
wavefunctions of Bogoliubov excitations of Bose conden- 
sate under the critical supercurrent in the presence of 
potential barrier with an arbitrary shape. With use of 
this result, we prove the partial transmission (i.e., the 
absence of perfect transmission ) of excitations in the low 
energy limit under the critical current. The low energy 
limit of wavefunctions of excitations does not coincide 
with that of the condensate wavefunction, as a result of 
the emergence of low energy density fluctuations local- 
ized around the potential barrier. 

Time-dependent Gross-Pitaevskii (GP) equation in 
dimensionless form is given by 

i^V(a:,t) = {-]^^ + U{x)^i;{x,t) + \^:{x,t)\^^{x,t). 

One can always recover the dimensional form by deflning 
a^dim = ix, tdim = ht/{gno), ?7dim(a;dim) = gnoU{x), and 
V'dimCa^dimjtdim) = y/uo'tp (x , t) , whcrc g is the strength 
of repulsive interaction, no is the condensate density far 
from the barrier, and ^ — h/y/mgriQ is the healing length. 
Setting the condensate wavefunction in the form of 

t) = e-''^*{*(a;) + [u(x)e-'^* - v*{x) e''*] } (1) 

and taking the terms up to first order with respect to 
u{x),v{x), we obtain the stationary GP equation 

L^{x) = 0, L = + U{x) \^{x)\' (2) 
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= e 



(3) 



for condensate wavefunction and Bogoliubov equations 

for wavefunctions of excitations. 

Wc assume that U{x) is short-ranged, and consider 
the solution which has the following asymptotic form: 



'u{x) 



u{x) 
—v{x)_ 



^{x — > ±00) = exp i^qx ^ ^ + const. j 



(4) 



The supercurrent q depends on the phase difference 
ip between two condensates separated by the potential 
barrier. The symbol "const." represents a trivial non- 
uniqueness of the phase factor. 

Setting ^'(x) = A{x)&^^^^ in eq. (2), one obtains 
2 

(5) 



dx \ dx 



= 0. 



(6) 



From (4), the boundary condition A(±oo) = 1 and the 
chemical potential ji = 1 + q^ /2 follow. From (4) and (6), 
we can obtain the expression for the phase of condensate 



<6{x) = q / = qx + q dx . .. „ 

and the expression for the phase difference 



1 



+00 



da; 



1 



A{xy 



1 



(7) 



(8) 



Here we fix Q{x = 0) to be zero. By eliminating d0/da; 
in (5), we obtain 



-1-1+^^. (9) 



Generally, eq. (9) has one or more solutions for given q, so 
if is more suitable than 5 as a parameter which specifies 
a state of the system uniquely. Therefore, we consider q 
and A as functions of ip, i.e., q{ip) and A(x, (p). q-ip curve 
corresponds to Josephson relation.*'' An example is 
available in Fig. 2 of ref. 13. The state where q takes an 
extremal value is the critical current state. 
Introducing 



S' = Me-"='+i;e"=', G 
eq. (3) can be rewritten as 

HS 



ue 



-ie 



ve 



ie 



ig_d_/G 
Adi\A 



eG, 



(10) 



(11) 



(12) 



S and G have the simple physical meanings; IV'P and 
ip/lipl of (1) can be written, respectively, as 



= A" 1 + - Re(Ge-'^*) 

L ^1 



(13) 



-i/it+ie 



1 + - Im(S'e-'<^*) 



(14) 



where higher-order terms of S and G are ignored. From 
these expressions, we can regard S/A and G/A as the 
phase and density fluctuations, respectively. It has been 
also pointed out by Fetter and Rokhsar.^^^ 

In the following, we consider the tunneling problem 
of excitations and hence we seek for the solution having 
the asymptotic form: 



S{x) 



{x — > —00) 
{x — > +00) 



(15) 



ki and ^2 arc real positive and negative roots of the 
dispersion relation e = qk + \yjk'^{k'^ + 4), and kx ~ 
e/{l + q) and ^2 — e/(— 1 + 9); respectively. It should 
be noted that |tp represents a transmission coefficient, 
while |fp is not a reflection coefficient. 

We construct the solution of the Bogoliubov equa- 
tions in the form of the power series in e : 



'S{x) 



n=0 



'S^"\x 
G(")(xi 



(16) 



Substituting it to the Bogoliubov equations, one obtains 
the recurrence relations: 



iq d fG^"" 



A dx 



(F + 2A2)G(")-^A('^) = 5(" 
Adx\ A J 



Q{n-1) 



(17) 



(18) 



From these equations, we see that ^("l G(") can be deter- 
mined from S^^\ G'*^-*, recursively. The general solution of 
eqs. (17) and (18) should be expressed in terms of four 
homogeneous and one particular solution(s): 



S'-^Xx) 
G(»)(x) 



G, 



= E 

j=I,II,III,IV 



'Sj{x) 



Gi%) 



(19) 



Here the last term in the right hand side does not exist 
if n = 0. 

First, we look for the general solution for n — 0, or 
equivalently, the homogeneous solutions {Sj,Gj). Since 
Bogoliubov equations always have the solution e = and 
(m, v) = (*, **),!«) one can take {Si, Gi) = {A, 0). So the 
solution of eq. (17) for n = becomes 

^G(o)dx 



S^°^ = CiA + CuAj^ ^^2\qAj^ 



A^ 



(20) 



by regarding G^"-* as an inhomogeneous term. Substitut- 
ing (20) to (18) for n = 0, we obtain the following equa- 
tion for G^: 

'h + 2A^ - gW = i^GiiA-^. (21) 
On the other hand, differentiating eq. (9) with respect to 
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ip, one obtains 



dA 
dip 



(22) 



Comparing (21) with (22), we find that A^ becomes a 
homogeneous solution of (21), only when the condition 
dg/d(/? = of critical current holds. Since A(x^±oo) = 
1, A^{x^±oo) = 0, i.e., A^{x) is localized near the po- 
tential barrier. As an example, see Fig. 1. Once we know 
this solution, other solutions can be obtained by meth- 
ods of reduction of order and variation of parameters. 
We thus obtain all homogeneous solutions as follows: 



Gill 



Sii 



-2\qAA-i 

Am 



-2igA, 



I 



-2iqA 



V 



Giydx^ 

dx' 



Here we have introduced the following function: 

rA^{x')dx' 



Aaix) 



r.'\3 



(23) 



/o M^' 

We note that the third solution (6*111, Gm) can be ex- 
pressed as {u,v) = {d^/dip,—d'^*/dip) by means of (7) 
and (10). We also note that (5n, Gu) and {Siy, Giv) are 
exponentially divergent solutions. 

Since we know all homogeneous solutions, the par- 
ticular solution (5^"\g|,"') in (19) can be obtained by a 
method of variation of parameters. It is 

dx' 



2iq 



dx" A^p{x") 



dx"'A{x"')G^-% 



A{x' 



S^^-^\x") 



(24) 



dx' 



dx"A{x")G'^''-^\x") 



- 2i9A(x) / 
Jo 



G^"V)da:' 
A{x'Y ■ 



(25) 



From now on, we assume that the barrier U{x) is 
even for simplicity. Accordingly, A{x) and A^{x) are also 
even. The asymptotic behavior of A^{x) is 



A^{x)^ l3e-^^^\^\ (a; ^±00), 



(26) 



where /3 is a certain constant. This asymptotics comes 
from the solution of (21) far from the barrier {U{x) ~ 
0,A{x) ~ 1). 

In the tunneling problem, wc arc interested in so- 
lutions free from exponential divergence. Accordingly, 
(^'("^G^"^) must be chosen to behave asymptotically as 
a polynomial of finite order in x. In general, {S'^\ Gp"^) 




Fig. 1. Localized density fluctuation solution A,p = dA/dip for 
a delta-functional barrier Uix) = Vq6{x) with Vq = 3.32. (We 
multiply by a constant in the figure.) 



itself diverges exponentially. Nevertheless we can con- 
struct a particular solution that docs not diverge expo- 
nentially by making a linear combination with (^n, Gn) 
or (5iv,Giv). 

Now, wc find (S'(i), G(i)) when a zero-energy solution 
free from divergence is chosen to be {S^°\ G(°)) = {Si, Gi) 
or (S'iii,Gni). For this purpose, we introduce the nota- 
tions Ai{x) := J^dx' A{x')A^{x'), ai :— ^i(+oo), :— 
^3(4-00), and ?7 := 0:1/0:3. Both Ax{x) and ^3(3;) are 
odd functions, and from eq. (26) they behave as 



Ai{x) sgna; aj - 



2^/T- 



(27) 



with i = 1,3. We note that 0:3 can be simplified by means 
of eq. (8): 03 = dx (^"^ - 1) = However, 

we hold this notation for a moment. 

When (5'("),G(°)) = {Si.Gi) = (.4.0), the solution 
that does not diverge exponentially is given by 



G«(x) 



s'^Xx) 
G'iXx) 



19 Jo Jo 



Sii{x) 
Gii(x) 

^G«dx\ 



1 
19 



A^ 



-2A, 



dx 



Ai - TjAs 
^1 



(28) 



+ ±C>C'/ 



— (x-|-7sgna;) ; ^ (x-|-i^sgnx) 

\q \ — q^ 



{\-n)l{2{\-e)) 

where 7 and v are certain constants. When (5*^"^ , G*^*^)) — 
(5*111, Giii) = (— 2iQ'AA3, A(^), on the other hand, we ob- 
tain 



-2A [ 
Jo 



p 
p 

^Aidx 
~A^ 



G^^\x) 



2iqA> ^i"'*''^ 



A3 



4iqAir, / dx 



/ 

Jo 



A1A3 - aias 



(29) 



"3 



-2j?(|x| + A) 



2g'(l + r?) 
1-02 



(|x| + «)' 



-(19(1 + ??))/(! - 9^)sgnx 



3/4 



J. Phys. Soc. Jpn. 



Letter 



where A and k are certain constants. Thus, non-divergent 
solutions up to first order in e are given by 



+ 0(e'), (30) 



with i = I, III, and their asymptotic forms at x 
become 

^2 



Sl°^''\x) 

Slf^''\x) 
-2iqa3 

V 

7= --7- 
1? 



1 + e 



sgn x + e 
i g(l - V) 

1 2 
1 — 



g2) 



a; + 7 sgn x 



ig(l - q^) 



\x\+\ 



A 



ri ig(l + 77) 



-A 



\q 1 - 

With the above solutions, we derive the transmission 
coefiicient t in (15). We expand (15) in 1 <C \x\ <C 1/fci, 
and expand coefficients by e such ss t = t^°^ + et^^^ + 
■ ■■ , f = f(°) + ef(^) H , and we obtain 



-K. 



±00 



(31) 



(32) 



(33) 



l + f(0) + e(f« + (^ + ^)x) 
iW+e(iW + ^x)+0(e2) 



(34) 



We then equate this expression with the asymptotic form 
of the general linear combination of solutions (31) and 
(32), that is. 



S{x) 



CiSl°'''\x) 



oioial 

r< '-'ill 



(x) 



-^iii" 



Expanding Cj and Cm as = C^°^ +eC^'' + --- , C, 



-2iqa3 

^(1) 



(35) 



(0) 

III 



(1) 
III 



III 



we can obtain simultaneous equations. 



and the solutions can be given by 



2qr] 



f r;) 



,2 ' 



<-"iii — 



/ - ri) 



(36) 
(37) 



It is obvious that < 1^^°^^ < 1 holds unless r/ = ig, 
thus the partial transmission in the low energy limit fol- 
lows. We note that S{x) must be expanded up to second 



order to obtain i^^^ and r 

Let us make sure that our expression reproduces the 

result in rcf. 6 for the delta- fimctional barrier U{x) = 
Vq5{x). The condensate wavefunction is given by^-' 
A{xf = 7(a;)2 + 7(3;) = ^/TV tanh(vT^(|a;| + 
a;o)). Since Aq/dip = at the critical point, q can 
be regarded as a constant for differentiation with re- 
spect to if. Accordingly, Ad A/ dip = ^d"f/dip = 7(1 — 
q^)dxo/d^ cos\i-\y/T^{\x\+xo)) 
sion, we can obtain rj — q^ + 7(0)^ = 
when Vb > 1, gc ^ 7(0) ^ so r; 

which is consistent with ref. 6. Though ref. 
6 gives the explicit expression only for the high barrier 
case, we have confirmed that our expression is exact irre- 



Using this expres- 
A(0)2. From ref. 6, 
Therefore 



spective of the height of barrier. We have also confirmed 
that r} = q occurs only when Vq = 0, i.e., no barrier ex- 
ists. Therefore, we expect that rj = ±q does not occur for 
a generic potential barrier, though we do not have the 
general proof for this criterion yet. 

Finally, we discuss the physical origin of the dis- 
appearance of anomalous tunneling under the critical 
current. RecaUing that {Si,Gi) and (5iii,Giii) corre- 
spond to (m,w) = (*,**) and (d'i' /dip, -d^* /dip), re- 
spectively, and using (35) and (37), we obtain 



lim 



oc 



g+r? dipy-"^*) ' ^^^^ 

i.e., lime_>o(w, w) 7^ i^,^*). This is in contrast to the 

case without supercurrent, in which the wavefunctions 
of excitations coincide with the condensate wavefunc- 
tion.^^^ In the absence of supercurrent, the diflferential 
equation for G^'^-' has the unphysical exponentially di- 
vergent solutions only, and G^°) cannot contribute to the 
wavefunction of excitations in the low energy limit. In the 
present case, on the other hand, the solution G\\\=A^p, 
which is localized near the potential barrier, contributes 
to the low energy wavefunction of excitations. Thus the 
presence of local density fluctuation near the barrier in 
the low energy limit is the origin of the absence of anoma- 
lous tunneling in the critical current state. 
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